LBNL-41309 



Spin- up of Solitary Pulsars: Signal of Phase Transition^ 



Norman K. Glendenning 

Nuclear Science Division & Center for Nuclear and Particle Astrophysics 
Lawrence Berkeley National Laboratory, Berkeley, California 

December 30, 1997 



G 
c 

X 
CD 
"O 
_C 

u> 
c 

CD 
CD 




log t (year) 



Neutron Stars and Pulsars 

November 17 to 20, 1997 
Rikkyo Univiersity, Tokyo, Japan 
To be published in Universal Academy Press 
(Organizers: N. Shibazaki and N. Kawai) 



1 This work was supported by the Director, Office of Energy Research, Office of High 
Energy and Nuclear Physics, Division of Nuclear Physics, of the U.S. Department of 
Energy under Contract DE-AC03-76SF00098. 



Spin-up of Solitary Pulsars: Signal of Phase Transition 



Norman K. Glendenning 

Nuclear Science Division & Center for Nuclear and Particle Astrophysics 
Lawrence Berkeley National Laboratory, Berkeley, California 

Abstract 

A phase transition in the nature of matter in the core of a neutron 
star, such as quark deconfmement or Bose condensation, can cause 
the spontaneous spin-up of a solitary millisecond pulsar. The spin-up 
epoch for our model lasts for 2 x 10 7 years or 1/50 of the spin-down 
time (Glendenning, Pei and Weber 1997). 

1 Introduction 

Neutron stars have a high enough interior density as to make phase tran- 
sitions in the nature of nuclear matter a distinct possibility. According to 
the QCD property of asymptotic freedom, the most plausible is the quark 
deconfmement transition. According to lattice QCD simulations, this phase 
transition is expected to occur in very hot (T ~ 200 MeV) or cold but dense 
matter. 

Since neutron stars are born with almost the highest density that they 
will have in their lifetime, being very little deformed by centrifugal forces, 
they will possess cores of the high density phase essentially from birth if the 
critical density falls in the range of neutron stars. However the global prop- 
erties such as mass or size of a slowly rotating neutron star are little effected 
by whether or not it has a more compressible phase in the core. In principle, 
cooling rates should depend on interior composition, but cooling calculations 
are beset by many uncertainties and competing assumptions about compo- 
sition can yield similar cooling rates depending on other assumptions about 
superconductivity and the cooling processes. Moreover, for those stars for 
which a rate has been measured, not a single mass is known. It is unlikely 
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that these measurements will yield conclusive evidence in the present state 
of uncertainty. 

Nevertheless, it may be possible to observe the phase transition in mil- 
lisecond pulsars by the easiest of measurements — the sign of Cl. The sign 
should be negative corresponding to loss of angular momentum by radiation. 
However, as we will show, a phase transition, either first or second order, 
that occurs near or at the limiting mass star, can cause spin-up during a 
substantial era compared to the spin-down time of millisecond pulsars. We 
sketch the conventional evolutionary history of millisecond pulsars with the 
addition of the supposition that the critical density for quark deconfinement 
falls in the density range spanned by neutron stars. 

As already remarked, with this supposition, the star has a quark core 
from birth but its properties are so little effected that this fact cannot be 
discerned in members of the canonical pulsar population. However, by some 
mechanism, usually assumed to be accretion from a companion during the 
radio silent epoch following its 10 7 year spin down as a canonical pulsar, 
the neutron star may be spun up. As it spins up, it becomes increasingly 
centrifugally deformed and its interior density falls. Consequently, the radius 
at which the critical phase transition density occurs moves toward the center 
of the star — quarks recombine to form hadrons. When accretion ceases, and 
if the neutron star has been spun up to a state in which the combination of 
reduced field strength and increased frequency turn the dipole radiation on 
again, the pulsar recommences spin-down as a visible millisecond pulsar. 

During spin-down the central density increases with decreasing centrifugal 
force. First at the center of the star, and then in an expanding region, the 
highly compressible quark matter will replace the less compressible nuclear 
matter. The quark core, weighed down by the overlaying layers of nuclear 
matter is compressed to high density, and the increased central concentration 
of mass acts on the overlaying nuclear matter, compressing it further. The 
resulting decrease in the moment of inertia causes the star to spin up to 
conserve angular momentum not carried off by radiation. The phenomenon is 
analogous to that of "backbending" predicted for rotating nuclei (Mottelson 
and Valatin 1960) and discovered in the 1970's (Johnson et.al. 1972, Stephens 
and Simon 1972) (see Fig. H) In nuclei, it was established that the change in 
phase is from a particle spin-aligned state at high nuclear angular momentum 
to a superfluid state at low angular momentum. The phenomenon is also 
analogous to an ice skater who commences a spin with arms outstretched. 



2 



Initially spin decreases because of friction and air resistance, but a period 
of spin-up is achieved by pulling the arms in. Friction then reestablishes 
spin-down. In all three examples, spin up is a consequence of a decrease 
in moment of inertia beyond what would occur in response to decreasing 
angular velocity. 

2 Calculation 

In our calculation, nuclear matter was described in a relativistically covariant 
theory (Garpman, Glendenning and Karant 1979, Glendenning 1985, Glen- 
denning and Moszkowski 1991) and quark matter in the MIT bag model 
(Chodos et.al. 1974). The phase transition occurs in a substance of two 
conserved quantities, electric charge and baryon number, and must be found 
in the way described in Ref. (Glendenning 1992). The moment of inertia 
must incorporate all effects described above — changes in composition of mat- 
ter, centrifugal stretching — and frame dragging, all within the framework of 
General Relativity. The expression derived by Hartle is inadequate because 
it neglects these effects (Hartle and Thorne 1968). Rather we must use the 
expression derived by us (Glendenning and Weber 1992a, 1992b). 

For fixed baryon number we solve General Relativity for a star rotating 
at a sequence of angular momenta corresponding to an equation of state that 
describes the deconfinement phase transition from charge neutral nuclear 
matter to quark matter. The equation of state is shown in Fig. |2| The 
moment of inertia as a function of angular momentum does not decrease 
monotonically as it would for a gravitating fluid of constant composition. 
Rather, as described above, the epoch over which an enlarging central region 
of the star enters the more compressible phase is marked by spin-up (Fig. |3|). 

3 Spin-up Era 

To estimate the duration of the spin-up, we solve the deceleration equation 
for the star with moment of inertia having the behavior shown in Fig. [3]. 
From the energy loss equation 
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Figure 1: Nuclear moment of in- 
ertia function of squared fre- 
quency for 158 Er, showing back- 
bending in the nuclear case. 



Figure 2: Equation of state for 
the fist order deconfinement phase 
transition described in the text. 



for magnetic dipole radiation we find 
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This expression reduces to the usual braking equation when the moment 
of inertia is held fixed. The braking index is not constant but varies with 
angular velocity and therefore time according to 



ftft 3/'ft + /"ft 2 
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where /' = dl/dfl and /" = dl 2 /dVt 2 . This holds in general even for a star 
whose internal composition does not change with angular velocity (inconceiv- 
able). In particular one can see that for very high frequency, the derivatives 
will be largest and the braking index for any millisecond pulsar near the 
Kepler frequency will be less than the dipole value of n = 3. 
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Figure 3: Moment of inertia 
of corresponding to a change of 
phase. Time flows from large to 
small I. 
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Figure 4: The time evolution of 
the braking index plotted over one 
decade that includes the epoch of 
the phase transition. 



The braking index is shown as a function of time in Fig. [| An anomalous 
value endures for 10 8 years corresponding to the slow spin down of the pulsar 
and the corresponding slow envelopment of a growing central region by the 
new phase. The two points that go to infinity correspond to the infinite 
derivatives of / at which according to (0), the deceleration vanishes. They 
mark the spin-up era. The actual spin-up lasts for 2 x 10 7 years or 1/50 of 
the spin-down time for this pulsar. This could be easily observed in a solitary 
pulsar and would likely signal a phase transition. 

4 Summary 

We have found the spin-up phenomenon in a pulsar model for which the 
change in phase occurs in a star near the maximum mass. The phase transi- 
tion can be first or second order as long as it is accompanied by a sufficient 
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softening of the equation of state. We have not quantified the meaning of 
"sufficiently soft". 

If no pulsar is observed to produce the signal, little is learned. Just as 
in nuclear collisions, failure to observe a signal does not inform us that the 
deconfined phase does not or cannot exist. 

This work was supported by the Director, Office of Energy Research, 
Office of High Energy and Nuclear Physics, Division of Nuclear Physics, of 
the U.S. Department of Energy under Contract DE-AC03-76SF00098. 
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